Addressing the complex phylogenetic relationship of the Gempylidae fishes using mitogenome data.

The Gempylidae (snake mackerels) family, belonging to the order Perciformes, consists of about 24 species described in 16 genera primarily distributed in tropical, subtropical, and temperate seas worldwide. Despite substantial research on this family utilizing morphological and molecular approaches, taxonomy categorization in this group has remained puzzling for decades prompting the need for further investigation into the underlying evolutionary history among the gempylids using molecular tools. In this study, we assembled eight complete novel mitochondrial genomes for five Gempylidae species (Neoepinnula minetomai, Neoepinnula orientalis, Rexea antefurcata, Rexea prometheoides, and Thyrsites atun) using Ion Torrent sequencing to supplement publicly available mitogenome data for gempylids. Using Bayesian inference and maximum-likelihood tree search methods, we investigated the evolutionary relationships of 17 Gempylidae species using mitogenome data. In addition, we estimated divergence times for extant gempylids. We identified two major clades that formed approximately 48.05 (35.89-52.04) million years ago: Gempylidae 1 (Thyrsites atun, Promethichthys prometheus, Nealotus tripes, Diplospinus multistriatus, Paradiplospinus antarcticus, Rexea antefurcata, Rexea nakamurai, Rexea prometheoides, Rexea solandri, Thyrsitoides marleyi, Gempylus serpens, and Nesiarchus nasutus) and Gempylidae 2 (Lepidocybium flavobrunneum, Ruvettus pretiosus, Neoepinnula minetomai, Neoepinnula orientalis, and Epinnula magistralis). The present study demonstrated the superior performance of complete mitogenome data compared with individual genes in phylogenetic reconstruction. By including T. atun individuals from different regions, we demonstrated the potential for the application of mitogenomes in species phylogeography.

The complete mitochondrial genome of the South African snoek Thyrsites atun (Euphrasén, 1791) (Perciformes, Gempylidae).

This is the first report on the complete mitochondrial genome (mitogenome) of South African Thyrsites atun (Euphrasén, 1791) and its phylogenetic placement within the Gempylidae family. The complete mitogenome of snoek is 16,494 bp in length and comprises 2 rRNAs, 13 protein-coding genes, 22 tRNAs, and one control region. Gene order is similar to that found in gempylids and other marine fishes. Reconstruction of Gempylidae phylogeny implies that the mitogenomes of snoek, black snoek Thyrsitoides marleyi, and snake mackerel, Gempylus serpens are closely related in evolutionary terms.

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip.

KEY MESSAGE: We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

Shallow seamounts represent speciation islands for circumglobal yellowtail Seriola lalandi.

Phenotypic plasticity in life-history traits in response to heterogeneous environments has been observed in a number of fishes. Conversely, genetic structure has recently been detected in even the most wide ranging pelagic teleost fish and shark species with massive dispersal potential, putting into question previous expectations of panmixia. Shallow oceanic seamounts are known aggregation sites for pelagic species, but their role in genetic structuring of widely distributed species remains poorly understood. The yellowtail kingfish (Seriola lalandi), a commercially valuable, circumglobal, epipelagic fish species occurs in two genetically distinct Southern Hemisphere populations (South Pacific and southern Africa) with low levels of gene-flow between the regions. Two shallow oceanic seamounts exist in the ocean basins around southern Africa; Vema and Walters Shoal in the Atlantic and Indian oceans, respectively. We analysed rare samples from these remote locations and from the South African continental shelf to assess genetic structure and population connectivity in S. lalandi and investigated life-history traits by comparing diet, age, growth and maturation among the three sites. The results suggest that yellowtail from South Africa and the two seamounts are genetically and phenotypically distinct. Rather than mere feeding oases, we postulate that these seamounts represent islands of breeding populations with site-specific adaptations.

The development of genome-wide single nucleotide polymorphisms in blue wildebeest using the DArTseq platform.

Blue wildebeest (Connochaetes taurinus taurinus) are economically important antelope that are widely utilised in the South African wildlife industry. However, very few genomic resources are available for blue wildebeest that can assist in breeding management and facilitate research. This study aimed to develop a set of genome-wide single nucleotide polymorphism (SNP) markers for blue wildebeest. The DArTseq genotyping platform, commonly used in polyploid plant species, was selected for SNP discovery. A limited number of published articles have described the use of the DArTseq platform in animals and, therefore, this study also provided a unique opportunity to assess the performance of the DArTseq platform in an animal species. A total of 20,563 SNPs, each located within a 69 bp sequence, were generated. The developed SNP markers had a high average scoring reproducibility (>99%) and a low percentage missing data (~9.21%) compared to other reduced representation sequencing approaches that have been used in animal studies. Furthermore, the number of candidate SNPs per nucleotide position decreased towards the 3' end of sequence reads, and the ratio of transitions (Ts) to transversions (Tv) remained similar for each read position. These observations indicate that there was no read position bias, such as the identification of false SNPs due to low sequencing quality, towards the tail-end of sequencing reads. The DArTseq platform was also successful in identifying a large number of informative SNPs with desirable polymorphism parameters such as a high minor allele frequency (MAF). The Bos taurus genome was used for the in silico mapping of the marker sequences and a total of 6020 (29.28%) sequences were successfully mapped against the bovine genome. The marker sequences mapped to all of the bovine chromosomes establishing the genome-wide distribution of the SNPs. Moreover, the high observed Ts:Tv ratio (2.84:1) indicate that the DArTseq platform targeted gene-rich regions of the blue wildebeest genome. Finally, functional annotation of the marker sequences revealed a wide range of different putative functions indicating that these SNP markers can be useful in functional gene studies. The DArTseq platform, therefore, represents a high-throughput, robust and cost-effective genotyping platform, which may find adoption in several other African antelope and animal species.

Population genetics of Southern Hemisphere tope shark (Galeorhinus galeus): Intercontinental divergence and constrained gene flow at different geographical scales.

The tope shark (Galeorhinus galeus Linnaeus, 1758) is a temperate, coastal hound shark found in the Atlantic and Indo-Pacific oceans. In this study, the population structure of Galeorhinus galeus was determined across the entire Southern Hemisphere, where the species is heavily targeted by commercial fisheries, as well as locally, along the South African coastline. Analysis was conducted on a total of 185 samples using 19 microsatellite markers and a 671 bp fragment of the NADH dehydrogenase subunit 2 (ND2) gene. Across the Southern Hemisphere, three geographically distinct clades were recovered, including one from South America (Argentina, Chile), one from Africa (all the South African collections) and an Australia-New Zealand clade. Nuclear data revealed significant population subdivisions (FST = 0.192 to 0.376, p<0.05) indicating limited gene flow for tope sharks across ocean basins. Marked population connectivity was however evident across the Indian Ocean based on Bayesian clustering analysis. More locally in South Africa, F-statistics and multivariate analysis supported moderate to high gene flow across the Atlantic/Indian Ocean boundary (FST = 0.035 to 0.044, p<0.05), with exception of samples from Struisbaai and Port Elizabeth which differed significantly from the rest. Discriminant and Bayesian clustering analysis indicated admixture in all sampling populations, decreasing from west to east, corroborating possible restriction to gene flow across regional oceanographic barriers. Mitochondrial sequence data recovered seven haplotypes (h = 0.216, π = 0.001) for South Africa, with one major haplotype shared by 87% of the individuals and at least one private haplotype for each sampling location except Port Elizabeth. As with many other coastal shark species with cosmopolitan distribution, this study confirms the lack of both historical dispersal and inter-oceanic gene flow while also implicating contemporary factors such as oceanic currents and thermal fronts to drive local genetic structure of G. galeus on a smaller spatial scale.

Identification and characterization of miRNAs transcriptome in the South African abalone, Haliotis midae.

Aquatic animal diseases are one of the most important limitations to the growth of aquaculture. miRNAs represent an important class of small ncRNAs able to modulate host immune and stress responses. In Mollusca, a large phylum of invertebrates, miRNAs have been identified in several species. The current preliminary study identified known miRNAs from the South African abalone, Haliotis midae. The economic and ecological importance of abalone makes this species a suitable model for studying and understanding stress response in marine gastropods. Furthermore, the identification of miRNA, represents an alternative and powerful tool to combat infectious disease.

Evaluation of de novo assembly technique in the South African abalone Haliotis midae transcriptome: A comparison from Illumina and 454 systems.

Next generation sequencing platforms have recently been used to rapidly characterize transcriptome sequences from a number of non-model organisms. The present study compares two of the most frequently used platforms, the Roche 454-pyrosequencing and the Illumina sequencing-by-synthesis (SBS), on the same RNA sample obtained from an intertidal gastropod mollusc species, Haliotis midae. All the sequencing reads were deposited in the Short Read Archive (SRA) database are retrievable under the accession number [SRR071314 (Illumina Genome Analyzer II)] and [SRR1737738, SRR1737737, SRR1737735, SRR1737734 (454 GS FLX)] in the SRA database of NCBI. Three transcriptomes, composed of either pure 454 or Illumina reads or a mixture of read types (Hybrid), were assembled using CLC Genomics Workbench software. Illumina assemblies performed the best de novo transcriptome characterization in terms of contig length, whereas the 454 assemblies tended to improve the complete assembly of gene transcripts. Both the Hybrid and Illumina assemblies produced longer contigs covering more of the transcriptome than 454 assemblies. However, the addition of 454 significantly increased the number of genes annotated.

Molecular systematics and biogeography of the circumglobally distributed genus Seriola (Pisces: Carangidae).

The genus Seriola includes several important commercially exploited species and has a disjunct distribution globally; yet phylogenetic relationships within this genus have not been thoroughly investigated. This study reports the first comprehensive molecular phylogeny for this genus based on mitochondrial (Cytb) and nuclear gene (RAG1 and Rhod) DNA sequence data for all extant Seriola species (nine species, n=27). All species were found to be monophyletic based on Maximum parsimony, Maximum likelihood and Bayesian inference. The closure of the Tethys Sea (12-20 MYA) coincides with the divergence of a clade containing ((S. fasciata and S. peruana), S. carpenteri) from the rest of the Seriola species, while the formation of the Isthmus of Panama (±3 MYA) played an important role in the divergence of S. fasciata and S. peruana. Furthermore, factors such as climate and water temperature fluctuations during the Pliocene played important roles during the divergence of the remaining Seriola species.

Genetic and Morphological Characterization of Freshwater Shrimps (Caridina africana Kingsley, 1882) Reveals the Presence of Alien Shrimps in the Cape Floristic Region, South Africa.

Morphological identification and molecular data (mtDNA COI) were used to resolve the taxonomic identity of a non-native freshwater shrimp in the Cape Floristic Region (CFR) of South Africa and to evaluate levels of genetic diversity and differentiation in the species' core natural distribution. The species was morphologically and genetically identified as Caridina africana Kingsley, 1882, whose main natural distribution is in the KwaZulu-Natal (KZN) Province, more than 1200 km from the point of new discovery. Subsequently, sequence data from natural populations occurring in seven rivers throughout KZN showed the presence of nuclear copies of the mtDNA COI gene (NUMTs) in 46 out of 140 individuals. Upon removal of sequences containing NUMTs, levels of genetic diversity were low in the alien population (possibly as a consequence of a bottleneck event), while varying levels of genetic diversity and differentiation were found in natural populations, indicating habitat heterogeneity, fragmentation and restricted gene flow between rivers. Following the present study, the alien shrimp has survived the Western Cape's winter and dispersed into a nearby tributary of the Eerste River System, hence posing an additional potential threat to endangered endemics. Understanding the biology of this alien species will aid detection and eradication procedures.

Domain repeats related to innate immunity in the South African abalone, Haliotis midae.

Molluscs predominately use the cellular defence system as the primary mechanism of defence against pathogenic infection, in which haemocytes play a pivotal role. Haliotis midae is a commercially important South African species that it is susceptible to bacterial pathogens, fungal and yeast infections in the farming environment. The current study aims to enrich the current knowledge regarding H. midae innate immunity by investigating the presence and evolution of domain repeats. The bioinformatics approach used in this study, detected five repeat families in the H. midae transcriptome. These repeats families include mixed alpha and beta (leucine-rich and ankyrin), spectrin repeats, beta-propellers (WD40) and alfa-structure repeat (TPR-like). The expansion of key gene families related to host defence may be important to abalone adaptation to life in a pathogen-rich environment.

Transcriptome profiles of wild and cultured South African abalone, Haliotis midae.

This report describes the use of pyrosequencing technologies to generate the first comparative analysis of de novo assembled transcriptome data from cultured and wild specimens of the South African abalone. The transcriptome data and database described here provide a significant genomic resource for abalone research. The data set annotated 11,240 genes, which matched genes with known functions in other species. A large number of transmembrane protein domains (4087) that may indicate a high portion of undiscovered gene receptors were identified. Further, we detected an interesting set of transcription factors (516) that are valuable candidates for participating in regulatory events in developmental (such as cell proliferation and differentiation) and reproductive processes.

Microsatellite cross-species amplification and utility in southern African elasmobranchs: A valuable resource for fisheries management and conservation.

BACKGROUND: Similarly to the rest of the world, southern Africa's diverse chondrichthyan fauna is currently experiencing high fishing pressures from direct and non-direct fisheries to satisfy market demands for shark products such as fins and meat. In this study, the development of microsatellite markers through cross-species amplification of primer sets previously developed for closely related species is reported as an alternative approach to de novo marker development. This included the design of four microsatellite multiplex assays and their cross-species utility in genetic diversity analysis of southern African elasmobranchs. As this study forms part of a larger project on the development of genetic resources for commercially important and endemic southern African species, Mustelus mustelus was used as a candidate species for testing these multiplex assays in down-stream applications. RESULTS: Thirty five microsatellite primer sets previously developed for five elasmobranch species were selected from literature for testing cross-species amplification in 16 elasmobranch species occurring in southern Africa. Cross-species amplification success rates ranged from 28.6%-71.4%. From the successfully amplified microsatellites, 22 loci were selected and evaluated for levels of polymorphism, and four multiplex assays comprising of the 22 microsatellites were successfully constructed, optimised and characterised in a panel of 87 Mustelus mustelus individuals. A total of 125 alleles were observed across all loci, with the number of alleles ranging from 3-12 alleles. Cross-species amplification of the four optimised multiplex assays was further tested on 11 commercially important and endemic southern African elasmobranch species. Percentage of polymorphism ranged from 31.8%-95.5% in these species with polymorphic information content decreasing exponentially with evolutionary distance from the source species. CONCLUSIONS: Cross-species amplification of the 35 microsatellites proved to be a time- and cost-effective approach to marker development in elasmobranchs and enabled the construction of four novel multiplex assays for characterising genetic diversity in a number of southern African elasmobranch species. This study successfully demonstrated the usefulness of these markers in down-stream applications such as genetic diversity assessment and species identification which could potentially aid in a more integrative, multidisciplinary approach to management and conservation of commercially important cosmopolitan and endemic elasmobranch species occurring in southern Africa.

Comparison of population genetic estimates amongst wild, F1 and F2 cultured abalone (Haliotis midae).

Haliotis midae is South Africa's most important aquaculture species. The reproduction cycle is currently not closed as many farms rely on wild-caught broodstock for seed production. However, there is an increasing interest in genetic improvement in commercial stocks, with a growing number of producers implementing selective breeding strategies. High throughput commercial production and mass spawning make it difficult to maintain breeding records; therefore, mostly mass selection is practised. The high fecundity and unequal parental contributions also often lead to increased levels of inbreeding. This study therefore aimed to assess the genetic effects of such breeding practices on commercial populations of H. midae. Using microsatellite loci, the genetic properties of a wild, an F1 and an F2 population were estimated and compared. Although there was no significant loss of genetic diversity amongst the cultured populations in comparison with the wild progenitor population, there was low-to-moderate genetic differentiation between populations. Relatedness amongst the F2 population was significant, and the rate of inbreeding was high. The effective population size for the F2 (±50) was also comparatively small with respect to the wild (∞) and F1 (±470) populations. These results suggest that farms need to give caution to breeding practices beyond the first (F1) generation and aim to increase effective population sizes and minimise inbreeding to ensure long-term genetic gain and productivity. This study also confirms the usefulness of population genetic analyses for commercial breeding and stock management in the absence of extensive pedigree records.

Transcriptome-wide single nucleotide polymorphisms (SNPs) for abalone (Haliotis midae): validation and application using GoldenGate medium-throughput genotyping assays.

Haliotis midae is one of the most valuable commercial abalone species in the world, but is highly vulnerable, due to exploitation, habitat destruction and predation. In order to preserve wild and cultured stocks, genetic management and improvement of the species has become crucial. Fundamental to this is the availability and employment of molecular markers, such as microsatellites and single nucleotide (SNPs). Transcriptome sequences generated through sequencing-by-synthesis technology were utilized for the in vitro and in silico identification of 505 putative SNPs from a total of 316 selected contigs. A subset of 234 SNPs were further validated and characterized in wild and cultured abalone using two Illumina GoldenGate genotyping assays. Combined with VeraCode technology, this genotyping platform yielded a 65%-69% conversion rate (percentage polymorphic markers) with a global genotyping success rate of 76%-85% and provided a viable means for validating SNP markers in a non-model species. The utility of 31 of the validated SNPs in population structure analysis was confirmed, while a large number of SNPs (174) were shown to be informative and are, thus, good candidates for linkage map construction. The non-synonymous SNPs (50) located in coding regions of genes that showed similarities with known proteins will also be useful for genetic applications, such as the marker-assisted selection of genes of relevance to abalone aquaculture.

Detection of molecular signatures of selection at microsatellite loci in the South African abalone (Haliotis midae) using a population genomic approach.

Identifying genomic regions that may be under selection is important for elucidating the genetic architecture of complex phenotypes underlying adaptation to heterogeneous environments. A population genomic approach, using a classical neutrality test and various Fst-outlier detection methods was employed to evaluate genome-wide polymorphism data in order to identify loci that may be candidates for selection amongst six populations (three cultured and three wild) of the South African abalone, Haliotis midae. Approximately 9% of the genome-wide microsatellite markers were putatively subject to directional selection, whilst 6-18% of the genome is thought to be influenced by balancing selection. Genetic diversity estimates for candidate loci under directional selection was significantly reduced in comparison to candidate neutral loci, whilst candidate balancing selection loci demonstrated significantly higher levels of genetic diversity (Kruskal-Wallis test, P<0.05). Pairwise Fst estimates based on candidate directional selection loci also demonstrated increased levels of differentiation between study populations. Various candidate loci under selection showed significant inter-chromosomal linkage disequilibrium, suggesting possible gene-networks underling adaptive phenotypes. Furthermore, several loci had significant hits to known genes when performing BLAST searches to NCBI's non-redundant databases, whilst others are known to be derived from expressed sequences even though homology to a known gene could not be established. A number of loci also demonstrated relatively high similarity to transposable elements. The association of these loci to functional and genomically active sequences could in part explain the observed signatures of selection.

Bioinformatic survey of Haliotis midae microsatellites reveals a non-random distribution of repeat motifs.

Recent studies have shown the non-random distribution of microsatellite motifs between genomic regions within a particular species. This study investigates such microsatellite distributions in the genome of the economically important abalone Haliotis midae, via a bioinformatic survey. In particular, the association of specific repeat motifs to coding regions and transposable elements is investigated. An understanding of microsatellite genomic distribution will facilitate more efficient use and development of this popular molecular marker. A bias toward di- and tetranucleotide repeats was found in the H. midae genome. CA microsatellite units were the most abundant repeat motif, but were notably underrepresented in genic regions where GAGT repeats predominate. Approximately 17.5% and 21% of the microsatellites showed gene and/or transposable element associations, respectively. This could explain the high genomic frequencies of particular motifs across the genome and may allude to a possible functional role. The data presented in this study are the first to demonstrate such non-random dispersal of microsatellites in abalone and support previous findings arguing in favor of non-random distribution of repeat motifs.

Differential growth-related gene expression in abalone (Haliotis midae).

The slow growth rate of Haliotis midae impedes the optimal commercial production of this most profitable South African aquaculture species. To date, no comprehensive effort has been made to identify genes associated with growth variation in farmed H. midae. The aim of this study was therefore to investigate growth variation in H. midae and to identify and quantify the expression of selected growth-related genes. Towards this aim, molecular methodologies and cell cultures were combined as a time-efficient and economical way of studying abalone transcriptomics and cell biology. Modern Illumina sequencing-by-synthesis technology and subsequent sequence annotation were used to elucidate differential gene expression between two sibling groups of abalone demonstrating significant growth variation. The expression of selected target genes involved in growth was subsequently analysed by quantitative real-time PCR (qPCR). Fast- and slow-growing abalone and in vitro primary haemocyte cultures treated with different growth-stimulating factors were used. The results obtained from transcriptome analysis and qPCR revealed significant differences in gene expression between large and small abalone, and between treated and untreated haemocyte cell cultures. Throughout in vivo and in vitro qPCR experiments, the up-regulation of genes involved in the insulin signalling pathway suggests that insulin may be involved in enhanced growth rate for various H. midae tissues.

Transcriptome characterization of the South African abalone Haliotis midae using sequencing-by-synthesis.

BACKGROUND: Worldwide, the genus Haliotis is represented by 56 extant species and several of these are commercially cultured. Among the six abalone species found in South Africa, Haliotis midae is the only aquacultured species. Despite its economic importance, genomic sequence resources for H. midae, and for abalone in general, are still scarce. Next generation sequencing technologies provide a fast and efficient tool to generate large sequence collections that can be used to characterize the transcriptome and identify expressed genes associated with economically important traits like growth and disease resistance. RESULTS: More than 25 million short reads generated by the Illumina Genome Analyzer were de novo assembled in 22,761 contigs with an average size of 260 bp. With a stringent E-value threshold of 10-10, 3,841 contigs (16.8%) had a BLAST homologous match against the Genbank non-redundant (NR) protein database. Most of these sequences were annotated using the gene ontology (GO) and eukaryotic orthologous groups of proteins (KOG) databases and assigned to various functional categories. According to annotation results, many gene families involved in immune response were identified. Thousands of simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) were detected. Setting stringent parameters to ensure a high probability of amplification, 420 primer pairs in 181 contigs containing SSR loci were designed. CONCLUSION: This data represents the most comprehensive genomic resource for the South African abalone H. midae to date. The amount of assembled sequences demonstrated the utility of the Illumina sequencing technology in the transcriptome characterization of a non-model species. It allowed the development of several markers and the identification of promising candidate genes for future studies on population and functional genomics in H. midae and in other abalone species.

Investigating the establishment of primary cell culture from different abalone (Haliotis midae) tissues.

The abalone, Haliotis midae, is the most valuable commodity in South African aquaculture. The increasing demand for marine shellfish has stimulated research on the biology and physiology of target species in order to improve knowledge on growth, nutritional requirements and pathogen identification. The slow growth rate and long generation time of abalone restrict efficient design of in vivo experiments. Therefore, in vitro systems present an attractive alternative for short term experimentation. The use of marine invertebrate cell cultures as a standardised and controlled system to study growth, endocrinology and disease contributes to the understanding of the biology of economically important molluscs. This paper investigates the suitability of two different H. midae tissues, larval and haemocyte, for establishing primary cell cultures. Cell cultures are assessed in terms of culture initiation, cell yield, longevity and susceptibility to contamination. Haliotis midae haemocytes are shown to be a more feasible tissue for primary cell culture as it could be maintained without contamination more readily than larval cell cultures. The usefulness of short term primary haemocyte cultures is demonstrated here with a growth factor trial. Haemocyte cultures can furthermore be used to relate phenotypic changes at the cellular level to changes in gene expression at the molecular level.